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We investigate the ground-state magnetic long-range order of quasi-one-dimensional quantum 
Heisenberg antiferromagnets for spin quantum numbers s = 1/2 and s = 1. We use the coupled 
cluster method to calculate the sublattice magnetization and its dependence on the inter-chain 
coupling J±. We find that for the unfrustrated spin- 1/2 system, an infinitesimal inter-chain coupling 
is sufficient to stabilize magnetic long-range order, in agreement with results obtained by other 
methods. For s = 1, we find that a finite inter-chain coupling is necessary to stabilize magnetic long- 
range order. Furthermore, we consider a quasi one-dimensional spin-1/2 system, where a frustrating 
next-nearest neighbor in-chain coupling is included. We find that for stronger frustration as well, a 
finite inter-chain coupling is necessary to have magnetic long-range order in the ground state, and 
that the strength of the inter-chain coupling necessary to establish magnetic long-range order is 
related to the size of the spin gap of the isolated chain. 

PACS numbers: PACS-key 75.10. Jm, 75.10.Pq, 75.40.Mg, 75.50.Ee, 



I. INTRODUCTION 

Low-dimensional quantum antiferromagnets have 
attracted much attention as model systems where strong 
quantum fluctuations may destroy magnetic long-range 
order in the ground state (GS).[l[ In particular, the 
one-dimensional (Id) quantum Heisenberg antifcrromag- 
net (HAFM) does not exhibit magnetic long-range or- 
der (LRO). In addition, there is a basic difference be- 
tween half-integer and integer antiferromagnetic Heisen- 
berg chains. @ While the Id HAFM with half-inte ger spin 
quantum number exhibits a gapless excitation spectrum 
and a power-law decay of spin-spin correlations, a faster 
exponential decay of spin-spin correlations - accompa- 
nied by a finite excitation gap A (spin gap) - is observed 
for integer-spin Id HAFM. However, it is known that for 
the Id spin-half HAFM a frustrated next-nearest neigh- 
bor exchange coupling may also open an excitation gap 
(for a more detailed discussion of Id spin systems, see 
[3]). For the formation of magnetic LRO in HAFM, 
the transition to two-dimensional (2d) lattices is cru- 
cial. The HAFM on 2d bipartite lattices exhibits mag- 
netic LRO at zero temperature and only competing in- 
teractions may destroy LRO (for a more detailed discus- 
sion of 2d spin systems, see 0, Q). In real materials 
we are often faced with the situation that the nearest- 
neighbor in-chain coupling J\ is dominant but an inter- 
chain coupling J± is also present. A very weak inter- 
chain coupling even seems to be a rare exception, see 
e.g. Rcf. @. Therefore, the study of quasi Id quan- 
tum HAFM's, i.e. systems where the in-chain couplings 
are larger than the inter-chain couplings, are - on one 
hand - of basic interest in connection with the dimen- 
sional crossover from one dimensions to two, and on the 
other hand, of interest for the interpretation of experi- 
ments. The quasi Id spin-1/2 HAFM has been studied 



in several papers @, H @, El El El El El El in recent 
years. A main focus of these studies has been on the 
estimation of the critical inter-chain coupling Jl where 
the transition between the phase with magnetic LRO and 
the magnetically disordered phase takes place. The an- 
swers given in the literature to this question are con- 
tradictory and not completely conclusive. While some 
papers find indications for a finite Jl d H EH j oth- 
ers find Jl = 0, E3, El El El E| which seems to 
be more plausible, given that the GS of the Id spin-1/2 
HAFM is not gapped. In particular, data obtaine d by 
the quantum Monte Carlo method (QMC) [III El Elf, 
which is precise for unfrustrated spin models, strongly 
support the result = 0. The behaviour of the quasi 
Id spin-1 HAFM might be different from its spin-1/2 
counterpart, but is less studied so far. Indeed, the ex- 
isting studies of the spin-1 case find indications for a fi- 
nite Jl @, El EE EE EE El which can be attributed 
to the spin gap between the singlet GS and the mag- 
netic excitations. In Ref. 0] a lower bound of Jl was 
estimated as Jl > 0.025Ji, whereas an upper bound 
Jl < 0.1892Ji was given in Ref. Recent QMC cal- 
culations [3 EE Elf yield Jl w 0.043 - 0.044 Ji, which 
is only about 10% of the Haldane gap A. 

To get experimental input for the theoretical work, 
materials with a quasi Id behavior and a small coupling 
ratio between inter- and in-chain exchange are needed. 
Experimen tally , materials such as Sr 2 CuQ3, 
Ca 2 CuO s USUI, Sr 2 V 3 9 M 1 ^BaCu 2 Si 2 7 [if, 
Sr 2 Cu{PC>4) 2 , Ba 2 Cu(PO^) 2 [23| arc quite good 
examples of quasi Id spin-1/2 antiferromagnets with 
"not-too-large" inter-chain coupling. The spin-1/2 anti- 
ferromagnet showing the most perfect Id behavior so far 
is SrCu{POi) 2 [||; this material has the smallest ratio 
fcsIjv/Ji ~ 6 x 10~ 4 between the Neel temperature Tjv 
and the in-chain coupling J\ . All of the above-mentioned 
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materials show magnetic LRO below Tjv, although their 
coupling ratio is remarkably small. Contrary to the spin- 
1/2 materials in some quasi-ld spin-1 antiferromagncts 
such as Ni(C 2 H s N 2 )2Nq2iCl0 4 ) @, H and other 



iV«-compounds, see Ref. [241 ]. no Neel LRO has been 



observed as measured down to very low temperatures. 

Another interesting system is the Id spin-1/2 HAFM 
with frustrated second neighbor interaction J 2 > 0. At 
J 2 = 0.2411Ji [2(| a transition to a dimerized state with 
a spin gap and an exponential decay of the spin-spin cor- 
relation occurs. Hence, for J 2 > 0.2411Ji the effect of 
the inter-chain coupling on the GS behavior may be dif- 
ferent from the case with J 2 < 0.2411J!. This problem 
has not been discussed in detail in the literature so far. 

In the present paper we apply the coupled cluster me- 
thod (CCM) [13, || to study the GS LRO of quasi Id 
HAFM with spin quantum numbers s = 1/2 and s = 1. 
This approach is a universal and powerful method of 
quantum many-body theory. Though the CCM is a fairly 
new method in the field of quantum spin systems, in re- 
cent years it has been developed to higher levels of ap- 
proximation which allows its application to quantum spin 
systems with more and more success (for recent reviews, 
see Refs. [29|,|3(|). In particular, the CCM has the advan- 
tage that it can be applied to frustrated quantum spin 
systems with arbitrary dimensions. Though, concern- 
ing the precision of the results, the CCM at the present 
level of approximation probably cannot compete with the 
QMC, it allows to find new results for frustrated systems, 
for which the QMC fails due to the so-called sign prob- 
lem. On the other hand, the comparison of CCM results 
with QMC data for the unfrustrated models can be con- 
sidered as a benchmark test of the CCM and is therefore 
of interest from an applied method point of view. 

The Hamiltonian of the quasi Id frustrated HAFM is 
written as 

H = {JlSi,n ' Si+i jT1 + J 2 Si^ n ■ Si + 2,ri) 

n i 

+ y] j±.s iin ■ Sj !n+1 . (i) 

i n 

The index n labels the chains and i the lattice sites within 
a chain n. The model is illustrated in Fig. Q] While for 
spin quantum number s = 1 we consider only the model 
without frustration (J 2 = 0), for s = 1/2, we will discuss 
both cases, J 2 = and J 2 > 0. 



II. THE COUPLED CLUSTER METHOD (CCM) 

In this section, the CCM formalism will first briefly 
be outlined. For further details the interested reader is 
referred to Refs. [H, H H H, H, . 
The starting point for the CCM calculation is the choice 
of a normalized reference or model state |<I>), together 
with a set of (mutually commuting) multi-configurational 
creation operators {C^} and the corresponding set of 



FIG. 1: Illustration of the quasi one-dimensional HAFM with 
the in-chain nearest-neighbor bonds Ji, the frustrating in- 
chain next-nearest neighbor bonds J2, and the inter-chain 
bonds J±, cf. the Hamiltonian JTJ. All bonds are antifer- 
romagnetic. 



their Hermitian adjoints {Cl}, 

(<S>\C+ =0 = C L |$) VL^O, C+^l (2) 
[C+, C+] = = [C L ,Cj}. (3) 

The operators (Cl) are defined over a complete set 
of many-body configurations denoted by the set-indices 
{L}. For the set {|$), C^} the CCM parametrization of 
the exact ket and bra GS eigenvectors l^) and (^1 of our 
many-body system are given by 



|$) , S = £ a L C\ 



(4) 



L#0 



(*| = {$\Se- s , S= 1 + J2<*lC l . (5) 

The CCM correlation operators, S and S, contain the 
correlation coefficients, and a^, which have to be cal- 
culated. Once these values are known, all the GS proper- 
ties of the many-body system can be derived from them. 
To find the GS correlation coefficients cll and cll, we 
simply require that the expectation value H = (^/\H\^) 
(GS energy) is a minimum with respect to the entire set 
{q,l, a^}, which leads to the GS CCM ket-state and bra- 
state equations 

($|C7e- 5 ffe s |$) = ; VL ^ (6) 
(<f>\Se~ s [H, C+]e 5 |$) =0 ; VL ^ 0. (7) 

For the spin systems considered herein, we choose the 
Neel state with spins aligned in the z-dircction as the ref- 
erence state. The reasoning behind this choice is evident 
for the unfrustrated case, [ [Hj]], since the Neel state is 
the classical GS for J 2 = 0. For the s = 1/2 case, the 
frustrated model (i.e. J 2 > 0) is considered below as 
well. Note that the classical GS is an incommensurate 
spiral state for J 2 > 0.25Ji. However, the quantum GS 
for J± = does not exhibit spiral ordering for values 
of J 2 less than J 2 ~ 0.5Ji, but it is rather a collincar 
state. [Ill [U, S3, EH Hence, the Neel state is an appro- 
priate reference state for < J 2 < 0.5 Ji as well. [3fJ] 
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To treat each side equivalcntly, we perform a rotation 
of the local axis of the up spins, such that all spins in the 
reference state align in the negative z-dircction. In this 
new set of local spin coordinates the reference state and 
the corresponding creation operators are given by 



|*} = IUU--->; C\ 



(8) 



where the indices i,j, k, ... denote arbitrary lattice sites. 
For the discussion of the GS Neel LRO, we have to calcu- 
late the order parameter (sublattice magnetization) M, 
which is given within the CCM scheme by 



M 



1 N 

i=l 



(9) 



The CCM formalism becomes exact if we take into ac- 
count all possible multispin configurations for the corre- 
lation operators S and S. However, in general, this is 
impossible to do in practice for a many-body quantum 
system. It is therefore necessary to use approximation 
schemes in order to truncate the expansions of S and S 
in Eqs. (|4I5[) in any practical calculation. A very general 
approximation scheme is the so-called SUB?i-to approx- 
imation. In this approximation, all correlations in the 
correlation operators S and S are taken into account, as 
long as they span a range of no more than m contiguous 
sites and contain only n or fewer spins. In most cases, 
however, the SUBn-n scheme is used (i.e., with n = to), 
and in these cases, it is referred to as the LSUBn scheme 
(for spin-1/2 systems). To find all the different (i.e. fun- 
damental) configurations entering S and S for a given 
level of SUB?i-n approximation, we use the lattice sym- 
metries. 

Although there is no theory available for how the re- 
sults of the SUBn-n approximations scale with n, there 
is nevertheless a great deal of experience in how to ex- 
trapolate the raw CCM SUBn-n data properly to n — > 
oo. [H [Ml, HI, IH, HI [H, H The best results for 
the extrapolation of the order parameter are obtained 
if the poor SUB2-2 data are omitted. Previously, for 
systems showing an order-disorder quantum phase tran- 
sition a leading 'power-law' extrapolation 
for the order parameter 



M(n) — ao + a\ 



(10) 



has been used successfully to determine the phase transi- 
tion points. In Eq. (fTU|) the leading exponent a-i is deter- 
mined directly from the SUBn-n data. Alternatively, as 
has been discussed recently in Ref. [39| , an extrapolation 
scheme with a fixed exponent but including an additional 
power in 1/n, i.e. 



M(n) = b 



1/2 



l 



(ii) 



can also be used to find the phase transition point. The 
extrapolation of the order parameter is illustrated in 
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FIG. 2: Extrapolation of the GS order parameter (sublat- 
tice magnetization) M scaled by the spin quantum number 
s for the pure unfrustrated s = 1/2 chain (i.e. J± = and 
J2 = 0). Fhe two methods of extrapolation corresponding 
to Eqs. (|10|l and (|lip are indicated by 'extrapolationl' and 
'extrapolation', respectively. 



Fig. [2] for one particular data set - namely, for the un- 
frustrated s — 1/2 chain with zero inter-chain coupling, 
i.e. at the expected critical point. Note that when using 
the extrapolation of Eq. (fit))) , the exponent a 2 for the 
considered data set is 02 = 0.414, which is a value not 
far from the fixed leading exponent 1/2 used in Eq. (|11[) . 
Below, we will use both extrapolation formulas to de- 
termine the critical inter-chain coupling J^_, which will 
naturally yield slightly different values for J£. We will 
use this difference as an estimation of the reliability of 
our results. 

Using parallel processing [35|, H(j we are able to use the 
CCM up to the SUB10-10 approximation for the spin- 
1/2 system, and up to the SUB8-8 approximation for the 
spin-1 system, which corresponds to the solution of more 
than 10 4 coupled nonlinear equations. 



III. RESULTS 

The results for the dependence of the order parameter 
M scaled by the spin quantum number s on the inter- 
chain coupling J± arc shown in Figs. [3H5J In those fig- 
ures, we show the raw CCM SUBn-n data used for the 
extrapolation, as well as the results for the extrapolation 
n — > 00, cf. Sect. IITI 

For the unfrustrated quasi Id HAFM with s = 1/2 we 
find that for J\ > J± > 0.2Ji, the variation of the or- 
der parameter M with J± is small, see Fig. [3] Only for 
Ji < 0.2Ji is the sublattice magnetization significantly 
diminished. As demonstrated in Fig. [31 it is evident that 
both extrapolation schemes lead to similar results. Using 
the extrapolation formula Eq. (|10[) . we obtain the criti- 
cal inter-chain coupling ~ 0.003 Ji at which the Neel 
LRO disappears . On the other hand, the extrapolation 
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FIG. 3: The dependence of the GS order parameter (sublat- 
tice magnetization) M scaled by the spin quantum number s 
on the inter-chain coupling J± for the unfrustrated s = 1/2 
quasi Id HAFM. The two methods of extrapolation corre- 
sponding to Eqs. pup and pip are indicated by 'extrapola- 
tion!' and 'extrapolation', respectively. 
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FIG. 5: The GS order parameter (sublattice magnetization) 
M scaled by the spin quantum number s in dependence on 
the inter-chain coupling J± for the frustrated s = 1/2 quasi 
Id HAFM. The two variants of extrapolation according to 
Eqs. pop and pip are indicated by 'extrapolationl' and 'ex- 
trapolation', respectively, a: J2 = 0.35Ji; b: J2 = 0.45Ji 
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FIG. 4: The dependence of the GS order parameter (sublat- 
tice magnetization) M scaled by the spin quantum number s 
on the inter-chain coupling J± for the unfrustrated s = 1 quasi 
Id HAFM. The two methods of extrapolation corresponding 
to Eqs. pop and pip are indicated by 'extrapolationl' and 
'extrapolation', respectively. 



based on Eq. (| 1 1 [1 leads to a finite but very small or- 
der parameter M Kt 0.01. These results are consistent 
with the conclusion that for the unfrustrated spin-1/2 
system, an infinitesimally small inter-chain coupling is 
sufficient to stabilize antiferromagnetic LRO. This state- 
ment supports the findings of Refs. @, UH, 0, d, Hi and 
are related to the gapless GS of the strictly Id s = 1/2 
HAFM. 

Next, we consider the unfrustrated s = 1 HAFM. Here 
the number of fundamental configurations in the CCM 
SUBn-n approximation is much larger than for s = 1/2, 
and we are able to present CCM data up to n = 8. Hence, 
the extrapolation for s = 1 is expected to be less reliable 



than for s = 1/2, for which CCM data up to n = 10 are 
available. Again, for J\ > J± > 0.2Ji, the variation of 
M with J± is small, see Fig. [U However, M/s is signif- 
icantly larger than as calculated for s = 1/2, indicating 
the decrease of quantum fluctuations in the square-lattice 
HAFM with increasing spin quantum number s. A strong 
reduction of the order parameter occurs for J± ~ 0.2Ji. 
Finally, M vanishes at a critical value J]_ > 0, i.e. a fi- 
nite inter-chain coupling is necessary to stabilize GS Neel 
LRO. However, both extrapolation schemes lead to dif- 
ferent numerical values for J^, namely J]_ ~ O.Ui using 
Eq. (|IDJ) but Jl - 0.02J 2 using Eq. JTTJ. Both values 
differ by approximately a factor of 2 from the QMC re- 
sults Jl « 0.043- 0.044Ji. This difference 
gives an indication of the accuracy of the CCM including 
SUBn-n data up to n = 8. 

Next, we consider the frustrated quasi Id s = 1/2 
HAFM for which QMC calculations are not possible. The 
classical GS for J2 > 0.25Ji is an incommensurate spi- 
ral state independent of the value of J±. In contrast to 
the classical model, for the strictly Id problem ( J± = 0) 
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FIG. 6: The critical inter-layer coupling JI in dependence 
on the frustration parameter J2 for the frustrated s = 1/2 
quasi Id HAFM. The two variants of extrapolation according 
to Eqs. (|10|l and are indicated by 'extrapolationl' and 
'extrapolation', respectively. 
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FIG. 7: The GS order parameter (sublattice magnetization) 
M scaled by the spin quantum number s for the frustrated 
s — 1/2 quasi Id HAFM in dependence on the frustration 
parameter J2 and for various strengths of the inter-chain cou- 
pling J ± . 



the quantum GS does not exhibit spiral correlations for 
0.25 Ji < J 2 < 0.5Ji 0, Ell, SI 51, but is gapped with 
a spin gap A strongly varying with J 2 . [42l. |44| However, 
in presence of an appreciable inter-chain coupling J±, 
spiral correlations may appear in the quantum model as 
well [13, SH|. Therefore, we restrict our calculations to 
J_l < 0.5Ji and to J 2 < 0.46Ji, where the CCM works 
well when based on a collincar reference state. The spin 
gap A of the Id problem (J± = 0) was found to be very 
small for J 2 < 0.4 Ji , but A increases rapidly between 

0.4 j! < j 2 < 0.5 j 2 . Hi m 

The results for M/s versus J± for J 2 = 0.35Ji and 
J 2 = 0.45 J 2 are shown in Fig. For J 2 = 0.35 Ji the 
critical inter-chain coupling is Jj_ ~ when using the ex- 
trapolation Eq. (fTU|) and is J^ ~ 0.014Ji when using the 



extrapolation Eq. (fTTj) . Knowing that for J? = 0.35Ji 
and J± = 0, the spin gap should be finite [26j but very 
small [U Ej], a zero or small is reasonable. On 
the other hand, for J 2 = 0. 45 Ji, where the spin gap 
for J ± = is A - 0.13Ji 0, [5| we already obtain 
quite a large value of J\ ~ 0.20Ji using the extrapola- 
tion of Eq. (fTfJ)) and JJ ~ 0.15Ji using the extrapolation 
of Eq. (JTTJ) ) . Interestingly, the ratio J^/A is seems to be 
larger than for the unfrustrated s = 1 HAFM. The varia- 
tion of JI with J 2 is shown in Fig.(S] Both extrapolation 
schemes yield qualitatively similar results. The variation 
of J£ with J 2 is quite similar to the variation of the spin 
gap with J 2 [42|, 0|. However, we observe a monotonic 
increase of the ratio R = J^/A from R ~ 0.8 to R ~ 1.6 
in the region 0.4Ji < J 2 < 0.46J 2 . 

The variation of the order parameter M/s with frustra- 
tion is illustrated in Fig. [71 where for the sake of clarity, 
only the data obtained using an extrapolation scheme 
Eq. pop are shown. As expected, frustration weakens 
the magnetic order and M becomes smaller with increas- 
ing J 2 . Obviously, for fixed but not too large inter-chain 
coupling Ji, a quantum phase transition between Neel 
LRO and a magnetically disordered phase can be driven 
by frustration. However, In similarity to recent findings 
for the quasi-2d Ji — J 2 model [38], it is likely that for 
stronger inter-chain coupling no magnetically disordered 
phase appears. 

Let us finally mention that the quantum phase tran- 
sition in the frustrated model is interesting from a more 
general point of view. For 0.2411Ji < J 2 < 0.5Ji, the 
model exhibits two ground state phases, each breaking 
different symmetries, namely (i) the rotationally invari- 
ant, spontaneously dimerized phase for zero (or small) Jj_ 
breaking the translational symmetry of the lattice, and 
(ii) the Neel phase, breaking the spin rotational sym- 
metry. A continuous transition between the dimerized 
phase and the Neel ordered phase is prohibited within 
the Landau theory. [46| , Hence three different scenarios 
are possible. First, that there is a (small) disordered fea- 
tureless spin-liquid phase between the dimerized and the 
Neel phase; second, that there is a first order transition 
between the dimerized and the Neel phase, and, third and 
most interesting, that the above-described transition is a 
candidate for a deconfined quantum critical point. (4fjj 
This question cannot be answered within the current ap- 
proach but deserves further consideration. 

To summarize, we find that the transition from the 
non-magnetic Id GS to the magnetically ordered 2d GS 
in quasi Id quantum HAFM's can be well described by 
the CCM if higher orders of approximation are used. Our 
results indicate that this transition, driven by the inter- 
chain coupling Jj_ , is related to the excitation gap of the 
strictly Id HAFM, i.e. at J± = 0. If the Id quantum GS 
is gapless, most likely the magnetic LRO sets in immedi- 
ately when J± is switched on, whereas for gapped GS's, 
a finite Ji is necessary to establish magnetic LRO in the 
GS. 
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